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A QUADRATIC PERFORMANCE INDEX FOR A VTOL 


AIRCRAFT PREFILTER MODEL REFERENCE 
ATTITUDE CONTROL SYSTEM 


Terrence D. Gossett* and Lloyd D. Corliss* 
Ames Research Center 


SUMMARY 


Several performance criteria methods for optimizing a prefilter model reference control 
system are reviewed, and the quadratic performance index is studied in detad. The prefilter model 
reference system consisted of a second-order model and a fourth-order follower. The form of this 
follower system was typical of that required for a VTOL aircraft in hover (e.g., the X— 14B). For 
the VTOL problem, the quadratic performance index was found to be the most versatile and 
comprehensive criterion. A system designer can iterate weighting matrices for various system states 
until the desired system response is obtained; that is, if the response is too oscillatory, the elements 
of the weighting matrices can be changed to produce a less oscillatory response. A major advantage 
of the quadratic performance index is that it can be minimized through the use of readily available 
digital computer programs. The solution of these programs provides a set of optimal gains for a 
given system configuration. Constraints on the use of the quadratic performance index are that the 
index must be expressed in terms of states and the control vector, all states must be observable, and 
only linear system formulations are acceptable. 

The shortcomings of applying the quadratic performance index are basically twofold: first, 
the process for selecting the weighting matrices is not well defined; and second, the parametric 
relationships of the weighting matrices to the resulting optimum gains are largely unknown. This 
report presents a simple graphical method for approximating the optimal gains. By this method 
asymptotes are generated for the inner-loop feedback gains in terms of the total outer-loop gain 
over a wide range of weighting. The asymptotes approximate the loop gains for the inner feedback 
and feedforward terms. 

The follower total outer- loop gain was found to be a useful reference parameter in 
constructing system asymptotes. The total outer- loop gain establishes the system “stiffness.” 
Analysis of the canonical form of the state equation and the Riccati equation for a single input 
single control system indicated that the follower total outer -loop gain could be analytically 
determined for type-zero or non-type -zero systems. However, only constant numerator, non-type- 
zero systems are studied in this report. The total outer-loop gain for these systems can be expressed 
in terms of the attitude weighting element of the Q matrix, the control weighting element of 
the R matrix, and the element from the control matrix, G. 


*U.S. Army Air MobilityR&D Laboratory, Moffett Field, California 


INTRODUCTION 


Historically, response feedback systems represent the conventional approach to augmentation 
of the inherent stabiUty of aircraft. However, the expansion of aircraft flight envelopes and 
variations in basic aircraft dynamic response due to changes in aircraft mass, inertia, and geometry 
have indicated a need for a more versatile class of flight control systems. The approach taken in this 
report is to consider a class of systems (i.e., model reference) that has inherent adaptability 
advantages and then to concentrate on the optimization of a fixed operating point of the system 
(the example used in this report is a VTOL aircraft in hover). The underlying intent is that one can 
then achieve adaptability by either discretely or continuously changing the system to comply with 
the chosen fixed operating points. The type of system felt to best fulfill the inherent adaptability 
requirement and also lend itself to the research investigation application is a model reference system 
of the pre filter model following type. This particular form has also been called the prefilter model, 
the model controlled system, the model reference adaptive system, or the model-in-the-system 
approach (refs. 1 to 5). 


The selection of an optimum model reference system for an aircraft application may be 
considered in two parts (as shown in fig. 1), both of which require the implementation of a 

performance criterion. The first portion of the 
problem involves the selection of a desirable 
model. This involves a man-machine interface 
for which the criterion of pilot opinion is most 
often used. This particular criterion is subjective 
but, nonetheless, seems well suited and, to date, 
the best available for such an interface. Several 
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Figure 1 .— Generalized model-follower representation. 


studies have been conducted utilizing pilot opinion to determine desirable hovering VTOL aircraft 
models in roll, pitch, and yaw (refs. 6 and 7). The roll model used in this study was one selected 
from these previous investigations (appendix A). 


The second portion of the problem, which is the focal point of this study, involves the 
optimization of the model-follower combination (fig. 1). For this link, which consists of a 
machine-machine interface, many performance criteria are available ranging from nonanalytic and 
rather subjective criteria, such as the comparison of time histories, to more analytical approaches, 
such as the quadratic performance index (QPI). Various criteria are discussed in appendix A. 
Nonanalytical criteria have the pitfall of not yielding a standard reference; that is, their “values” 
may vary from observer to observer. To develop a mean standard reference for such criteria requires 
the use of statistical methods. While a nonanalytical criterion, such as a pilot opinion, seems 
appropriate for the pilot-model link, it is deemed that an analytical criterion can and should be used 
for the model-follower link. This study, then, is concerned with the latter consideration, that is, the 
determination and interpretation of a satisfactory analytical performance criterion for the 
model-follower optimization. 


Having specified a prefilter model reference system through the optimization of these two 
links, the final test for the acceptability of the system rests with the pilot. Hence, the optimization 
of the follower through analytical means should not unduly influence the pilot’s appraisal of the 
system. In short, the pilots assessment of the overall system should hopefully be consistent with his 
assessment of the model alone. 
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For the model-follower optimization, several possible performance criteria were reviewed 
(refs. 8-13). These criteria included Integral Absolute Error, Integral Square Error, Integral Time 
Absolute Error, Butterworth, and the Quadratic Performance Index (QPI). The QPI, which may be 
minimized by solving the matrix Riccati equation, was found to be the most desirable criterion 
(refs. 5, 6, and 9). (See appendix A.) While reference 5 established a good basis for the application 
of the QPI, certain shortcomings of this technique were cited: (I) The process for selecting the QPI 
weighting matrices is not well defined. (2) Parametric relationships of these wei^ting matrices to 
the optimum gains produced are largely unknown. 

Alleviation of the latter limitation would enhance the use of QPI in practice on a synthesis 
problem. If parametric relationships are known, the designer could make an “on the spot” 
alteration of the system gains in an optimal manner. Presently, to achieve such an optimal alteration 
of gains for a given system requires a complete iteration of the Riccati equation solution. 
EstabUshment of a set of parametric relationships would in turn provide further insight on the first 
shortcoming cited above, that is, the selection of the weighting matrices. Hence, for a specific class 
of system, the main objective of this report is to relate analytically and graphically the effects of 
variations in the QPI weighting elements on the optimum gains they produce. While the 
relationships developed in this report are with reference to the X— 14B VTOL model reference 
system, the technique need in no way be restricted to such a system. 


GUIDELINES FOR INVESTIGATION 


For the purposes of this study it will be assumed that the X— 14B system will consist of a 
three-axis (roll, pitch, and yaw) prefilter model reference system that will utilize a general purpose 



Follower 


Figure 2.— Simplified roll axis model-follower diagram. 


digital computer and an autopilot for 
which the roll axis is shown on figure 2. 
The guidelines for investigation of 
performance criteria that follow were 
definitely influenced by this assumption. 
The general guidelines and procedures 
that will be observed in this study of 
performance criteria are outlined below: 

The form of the linearized lateral 
equations of motion characterizing 
hovering vehicles that will be used in this 
report is taken from reference 1 4. 






(■• ») 
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For hover, Uq = 0 and also for the simphfying, but generally valid, assumption that 
Ixz = Np/Iz = Yf/m = I^/Ix “ Ly/^x “ these equations reduce to: 



This report uses a single-axis roll transfer function for hover; hence it is necessary to find the 
transfer function for Therefore, Y^Jm, Lg^/m, and will b® considered to be 

accelerations from roll input only. In hover, the roll control generally causes only small side forces 
or yawing moments; hence Y^Jm and will be set to zero. Then from the above equation, 

0^/83 is found to be: 



The roll acceleration due to reaction nozzle input, Ls^/Ix, on the modified X— 14 is estimated to be 

4.585°/sec^/deg, and the measured value for Lp/lx for the X— 14 was — 0.45sec‘. Thus the 
transfer function that will be used for <^»j\/8a for hover is 

_ 4.585 

6^ s(s + 0.45) 

In hover, the X— 14 obtains its roll moment from reaction nozzles that are driven by a 10-cps 

actuator (^3 = 0.5) with a dc gain of 1 .0°/deg. Hence the actuator transfer function is 

<5 2 

^ ^ (62.8)^ 

^ s2 + 2(0.5) (62. 8)s + (62.8)^ 

Therefore, the plant representation of the basic aircraft is 

^ ^ (62.8)^(4.585) 

[s2 + 2(0.5) (62. 8)s + (62.8)2 ]s(s+ 0.45) 

^ 18,100 

s(s3 + 63.3 s 2 + 3978s + 1775) 
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Three other guidelines were established. Only the prefilter model reference system will be 
considered. The model in the examples will be second order with values of com = 2 rad/sec 
and “ 0-^ (which corresponds to a satisfactory hover handling qualities model as reported in 
ref. 7). The error to be minimized will be the difference between the corresponding responses of the 
model and the follower (e.g., the difference between the model and the follower attitudes; see 
fig. 2). 


ANALYTICAL INVESTIGATIONS 


Simply stated, the model following concept requires that the response of the follower match 
the response of the model, that is, in figure 2, a matching of to 0m- This implies then that any 
difference in the responses of the follower and the model constitutes an error signal. It is this error 
signal and its minimization that forms the basis of the various performance criteria discussed herein. 

One of the major objectives of this research was to better understand the effect of 
performance criteria on the system optimums that they produce. The search for a proper criterion 
was evolutionary in that the simpler criteria were explored first, deficiencies and merits were noted, 
and then more complex criteria were considered. 


QPI Criteria 

The quadratic performance index (QPI) was felt to be the most comprehensive and 
satisfactory criterion (refs. 5, 6, 9, and 15; see appendix A.) It can be minimized by solving the 
matrix Riccati equation with the aid of a digital computer and the programming package Automatic 
Synthesis Program (ASP) (ref. 9) or the more recently developed Fortran IV version of Automatic 
Synthesis Program (FASP) (ref. 16). The expressions used in this system optimization are given 
below; the terms are defined in the Symbols. These relationships are explained and developed in 
appendix B. 


V = / (XTrTqhx + uTRU)dt 

Quadratic performance index 

(1) 

P = 0 = PGR-* gTp - PF - fTp - hTqh 

Steady-state Riccati equation 

(2) 

U* = - R-‘ gTpx = K*X 

Qptimum control 

(3) 

Vmin =/xT(hTqh + PGR‘‘ gTP)X dt 

Minimized QPI 

(4) 

X = FX + GU 

System equation 

(5) 


The system given by equation (5) may be optimized, that is, the QPI may be minimized, 
through the solution of the Riccati equation. The matrix P is the solution to the matrix Riccati 
equation and is used in generating the optimal control U*. 
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It was recognized that various aspects of the application of the QPI are not easily understood; 
hence much of this investigation was directed toward the problem of gaining a better understanding 
of the optimum gains produced by this technique. This problem centers primarily on understanding 
the import of the weighting matrices Q and R and on developing a way of relating these Riccati 
optimums to other existing techniques. It is felt that the most effective way to explain the 
application of the QPI criterion is to present an example. 

Example for QPI 

This section will present the mathematical formulation of the optimal control problem 
described by equations (1) through (5) for a symbolic example. Then an analytical expression will 
be developed through the use of the matrix Riccati equation for one of the optimal gain terms. 

The example for this study is a prefilter model following system with the second-order model 
and fourth-order follower of the form shown in figure 3. Figure 3 is a generalized form of the 
system to be used on the X— 14B. The QPI weighting will be on the cost of control and on either 
the matching of attitudes or of attitudes and rates. The state equation for this example and the 
definition of the various matrices used in the QPI optimization procedure are given below. Further 
detail of the state space modeling and solution techniques appears in appendix B. 



Figure 3.— Detailed roll axis model-follower diagram. 


The nonhomogeneous state equation, that is, 

X = FX-i-GU (5) 
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when expanded for the generalized fourth-order follower system of figure 3, yields 


d 

dt 
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(5a) 


where 5 is a step input. For ease of illustration, equation (5a) may also be written in the following 
partitioned form: 



‘[Xf]' 


'[Ff] 

[0] 

[of 


"[Xff 


~[G{] 

d 

dt 

[Xml 


[0] 

[M] 

[B] 


[Xml 

+ 

[0] 
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[0] 
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(5b) 


Relating equations (5a) and (5b) yields the following subsystem matrices: 
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Fourth-order unaugmented follower matrix 


Second-order model matrix (step input) 
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0 
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Control matrix 


[Gf] 

[ 0 ] 

[ 0 ] 


§41 

0 

0 
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In addition, the matrices of the QPI expression (eq. (1)) for this example are given below. For 
convenience an error vector Y is defined as: 

Y = HX Error equation 

The matrix H is selected so as to yield as the error the combination of the states to be minimized 
by the QPI. For example, if the error Y is taken to be the differences between the model and 
follower attitudes and rates, an H of the following form is required: 


which yields 


H = 


-1 


0 0 0 1 


0-100 


0 0 


1 0 J 


Y 


(♦m - V 


The weighting matrices of the QPI integrand, which are also selected by the user, are in this study 
selected as shown below. 


Q = 



0 n 




R = r (a scalar) 



Figure 4.— Matrix representation of roll axis model-follower. 


Figure 4 is a structured diagram of 
the system as optimized by the QPI 
(eqs. (l)-(5)) and involving the above 
matrices; U* is the optimized form of the 
control U. Relating this optimized form to 
the example in figure 3 shows that the 
control U may be written as 


U = KX = -ko0A - ki 0A - ^2 '^ A ~ + ^5^ 

which, when optimized as given by equation (3), becomes 
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u* 



42 




P44 P45 P46 P47] 
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^A 

•I'm 
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(3a) 


where pjj is the ij element of the P matrix. Comparison of the above two expressions for u and 
u* reveals that in the optimized case ko = (g 4 i/r)p 4 i, kj = (g 4 i/r)p 4 2 , etc. Thus solving for the 
appropriate p elements provides the optimal K gains, which are defined as K*. 

Before proceeding with the solution of these p elements, one additional system expression 
will be constructed, namely, the generalized optimum follower transfer function 0A(s)/^m(s)- This 
expression will illustrate the similarities between the p elements, the K* gains, and the follower 
transfer function coefficients. 

This transfer function can be derived by substituting equation (3a) into line 4 of 
equation (5a), which yields 


_ N(s) 
(|>j„(s) " D(s) 


-(g4lVr) p^^s 

+ [(g4i^/r)p^^ - 

^ [(g41^/r)P42 - ^42^S 


+ P47[5(s)/<J>^(s)]f 

^ + [(g4iVr)p^3 - f43ls2'^ 

[(g4i^r)p^j - y 


( 6 ) 


The numerator of equation (6) can be further expanded by utilizing line 6 of equation (5a), that is. 
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which finally yields a numerator of the form 
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the generalized fourth-order optimum follower transfer function becomes 

4,^(s) NgS^ + N 2 S + 

+ DgS^ + D2S + Dj 


( 8 ) 


This form will be used in the study ; equations (7a) through (7g) are the key expressions showing 
the relationship between the p elements, the optimal gains K*, and the optimum follower transfer 
function coefficients. 


An expression for the optimum gains will now be determined in terms of the system and 
weighting matrices. The steady state Riccati equation (eq. (2)) can be constructed term by term; 
that is, 


0 = PGR-‘ gTr - PF- pTp - hTqh 


(2a) 


where 
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Since P is symmetric, the following is true; 
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In the Riccati equation, the elements of the matrices G and F are a function of the system 
under investigation and hence are known, and the elements of the matrices Q, R, and H are 
assigned by the designer. Thus only the elements of the P matrix (e.g., a 7 X 7 in the present 
example) must be solved for. However, it is not necessary to evaluate all 49 p elements because 
the P matrix is symmetric; this leaves 28 distinct p elements. Although, given a specific system, 
the computer conveniently provides a numerical solution for the p elements, information 
regarding the parametric relationships between the elements of the matrix P and the elements of 
matrices G, F, Q, R, and H is lost. Thus pursuing an algebraic approach, inspection of 
equation (3a) reveals that only the fourth row p elements need be evaluated to optimize the 
system. Accordingly, a generalized expression involving each of the fourth row p elements can be 
extracted from the Riccati equation, which is a function of G, F, Q, R, H, and also some of the 
other p elements. For example, the 2—2 element of equation (2a) yields an expression 
containing p 4 j: 


(-g4 1^/r)P4 2^ +2 (Pi 2 +f42P42) + Q22 =0 

This particular expression is a function also of the element Pi 2 and thus proves not to be a useful 
expression for solving for P 42 . In fact, only one element, P 41 , can be conveniently evaluated by 
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this algebraic method. Nonetheless, the expression for P 41 is a parametric one in terms of the 
system and the weighting matrices, thus satisfying one of the objectives of this report. The 
element p 4 1 also leads to the determination of one of the more “significant” optimal gain terms, 
as will be subsequently shown. Solving for P 41 is easily achieved through evaluation of the 
one-one element of the steady-state matrix Riccati equation, (2a); that is. 


(—§ 4 1 ^ /r)P 4 1 ^ 2 f 4 1 P4 1 Qi 1 “ 0 


Solving for P 4 1 yields 
r£ 


41 
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11^41 
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( use - for 
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f4i < 0\ 

^41 > 0/ 


Substituting this value for P 4 1 into equation (3a) yields ko , which is the outer-loop augmenting 
feedback gain term that operates on 0 ^ (see fig. 3) : 


-1 “+1 
^ SuiPl 


’ 41*^41 


>41 


+ 


+ 


qilg41^ \ 

^^41^ / 


Thus the optimum outer-loop augmenting feedback term for this example is a sole function of qj , , 
r, f 4 1 , and g 4 1 . The follower system total outer-loop gain term (i.e., the transfer function 
coefficient Dj ) is 


1 /2- 


*1 ^41^0 "'41 


- f,., = -f, , + f 


41 


■41 


1 ± 




1*41 


(9) 


rf 


41 


In non-type-zero open-loop follower systems (i.e., a system with one or more open-loop poles at the 
origin) this expression can be simplified. A non-type-zero system results if f 4 1 = 0 in figure 3 

ko = (Qi i/r)‘/^ 

and equation ( 9 ) becomes 

Di =g4i(qii/r)‘/^ (10) 

This relationship indicates that for a non-type-zero fourth-order system, the optimum outer-loop 
gain term is a function only of g 4 1 , qi i , and r. 


N Order Example of QPI 

Equations (9) and (10) were derived for a fourth-order follower system and P 41 was 
determined exactly. Before the other K* gain terms are discussed, a more generalized expression 
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Figure 5.— N-o-der type-zero follower system. 


of ko involving will be presented. 
Consider the n order type -zero system 
shown in figure 5. Through a development 
similar to that used for the fourth-order 
case, the optimum gain ko and the total 
outer-loop gain Dj can be determined. 
With the follower system 
matrix Ff represented in its canonical 
form, the element of F which enters into 


the optimized ko expression is fm- Correspondingly, the element of G is g^i , the element of Q is 
q, , , and the element of R is the scalar r. The ko expression can then be written as 


kn = r p 
“ ^nl^ni 


ni 


g. 


ni 


1 ± 1 + 


^11 ^ni' 


1 / 2 - 


r ^ 2 

ni 


Similarly, the total outer-loop gain term Dj is 

1 ^ni 0 ni ni ni 


1 ± 1 + 


fm ^ / J 


(9a) 


where 


n 


f = FB. 
ni • , 1 
1=1 


Again, for a system with fm =0 (non-type^/ero), equation (9a) reduces to 

= gni(9i 1 

Table 1 summarizes the analytical expressions for the n order case. 

TABLE 1.- ANALYTICAL SUMMARY - N ORDER OPTIMAL SYSTEMS 


Outer-loop 
augmenting 
gain ko 


Zero-Type Systems* 


ni 

T 

"ni 


1 ± 1 + 




1/2 


rf 


ni 


Non-Type-Zero Systems 
1/2 



Total outer- 
loop gain 
(TOLG = D, ) 


-f + £ 
ni ni 


1 ± 1 + 


^ll^nl' 


rf 


nl 


1/2 



* Use — for f 4 j less than 0. 

Use + for f 4 1 greater than or equal to zero. 
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Since the X-14 airplane equations are non-type-zero, fourth-order non-type-zero systems are 
considered in the remainder of this report. Hence, the outer- loop augmenting gains will be 
(qi 1 /r)‘ , and the total outer-loop gains will be , (qj i /r)‘ 

As indicated previously, unlike P 41 the outer p elements, and thus the other model 
feedforward and follower feedback terms, are not as easily expressed analytically through the use of 
the Riccati equation. Thus a graphical method was developed that utilizes the optimum outer-loop 
gain, which is the transfer function coefficient D, , and relates all the other optimum transfer 
function coefficients to it. 

The purpose of this graphical investigation is to show the relationship between Di (hence, fni, 
8 ni > <li 1 > r) remaining system optimum loop gains. Only the fourth-order follower 

system with constant numerator will be treated here although this development does not preclude 
its use for other order systems. Initially for illustrative purposes a simple system whose transfer 
function is 1/s^ will be considered for both the QPI attitude weighting case and then for the 
attitude and rate weighting case. Through optimization of this system for several Q and R matrices, 
graphical asymptotes can be constructed which, under certain conditions, become the asymptotes 
for other fourth-order systems. Further, for the attitude weighting only case, these asymptotes 
represent the Butterworth filter form coefficients. After the parametric asymptotes are established 
for the 1/s^ case, the graphical technique will be used to optimize the X— 14B system. 


Attitude Weighting Only 

To show the relationship between g 4 ,, q,i, r, and the inner-loop terms, first consider the 
preceding generalized fourth-order follower system in figure 3 for the specific case 
f4i~f42~f43~f44~0 (i-C-, transfer function of l/s**). The example with this follower does not have 
much practical significance, but is used to show that the following graphical technique obtains the 
Butterworth functional relationships that are known to exist for this case when the QPI is 
optimized (refs. 5 and 17). 

The optimal gains for the l/s"* follower can be determined parametrically without the aid of a 
computer by utilizing a direct solution detailed in reference 15. This method centers on the 
algebraic solution of the expression 


Is - (F - GK) 


-Is - (F - GK)^ 


-1 T 

Is - F GR G 

T T 

-H QH -Is - F 


After the determinants in the above equation are expanded and coefficients of like powers of s are 
equated from the left and right side of the equation, the following optimal gains result; 

ko = (Qi 1 /r)‘ 

ki =(2.6/g4//^)ko3/'‘ 
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^2 ~ (3.4/g4 1 * )ko * 


ka =(2.6/g4i"/^)ko>/^ 

For the I/s'* example, the optimal transfer function of the form of equation ( 8 ) simplifies to 


V(s) 


N 3 S + N 2 S + Nj 


d) (s) 4 , 3 . , 2 , , 

m s + g^kjS + g 4 ik 2 S + git^k^s + g^^kp 

Hence the above optimal gain expressions forko, kj, k 2 , and ka yield the following transfer 
coefficients: 

Di = g4 1 ko 


D 2 =g4iki =2.6 D,3/4 
Da =g4ik2 =3.4 D,‘/2 
D 4 =g 4 ,k 3 =2.6 D,'/'* 


Note that for this case the exponents and numerical coefficients are Butterworth. 

Although the analytical technique yields the parametric relationships for the 1/s^ follower, 
it may be tedious to apply in other cases. For those cases, a simple-to-apply graphical technique has 
been developed. The technique is based on the observation that the parametric relationships as 
calculated analytically for the 1 /s^ plant are of the general form 

Dn = bD,"^ 

This equation, when represented in log-log coordinates, is represented by a straight line, where 
constants b and m represent the intercept and slope of the line. 

Now pursuing the graphical approach, the I/s'* follower system was optimized through the 
solution of the Riccati equation for several values of qj 1 and r. The p elements for the special case 
of qi 1 = 10 , qj 2 = 0 , r = 0 . 01 , and 84 1 = 1 (i.e., qi 1 /r = 1000 ) resulted in the following optimum 
control (see eq. (3a)): 

u* = -100[0.316 0.348 0.192 0.062 -0.069 -0.042 -0.092JX 

where, from this computer solution, the total outer-loop gain term, Di , is 

Dj =g4iko = 31.6 

This can also be obtained analytically by utilizing equation ( 6 ); that is, 

r P 41 “ 841 y r j 


15 




Figure 6.— Optimum denominator terms for 1/s^ plant 
(attitude weighting only). 


A plot of this computer solution plus 
solutions for other Dj terms yielded the 
log-log plot (fig. 6 ) of the individual 
denominator terms D 4 , D 3 and D 2 plotted 
versus Dj . The values of g 4 , (q, , /r)' 
plotted were 3. 1 62, 3 1 .62, 3 1 6.2, and 3 1 62 
with g 4 1 = 1 and g 4 1 =10. This log-log plot 
might well have beenp 4 4 , P 43 , and P 42 
plotted versus P 4 , and the same 
straight-line functional relationship would 
still exist. However, the quantities used in 
figure 6 are preferred for substitution into 
the optimum transfer function of 
equation ( 8 ); that is, the optimum transfer 
function coefficients for any qi 1 /r can be 
determined from this representation. 


Figure 6 indicates that the relationships between the optimum total outer-loop gain and all of 
the optimum inner-loop gains can be exactly characterized by the equation shown earlier. 


For example. 


Dn = bD,"i 


D 2 = 2.6Di°-^5 


which, as anticipated, is the same parametric expression previously generated by the direct solution. 



I 10 100 1000 


D| =(g|/f)P4rq4|(P|/r)'^2 

Figure 7.— Optimum numerator terms 
for I/s'’ plant (attitude weighting 
only). 


In figure 7 the individual numerator terms in equation (8) 
(Ni, N2, and N3) are plotted against the same denominator 
term Dj . This log-log plot indicates the same straight-line 
functional relationship seen in figure 6. 

Six functional relationships for this system evolve from 
figures 6 and 7 : 


D 2 =2.6D,°-” 

D3 = 3.4Di°-5° 
D4 = 2.6Di“-2s 
Ni = Di 

N 2 =0.75Di“-^« 
N 3 = 0.289D,°-^“ 
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These significant relationships allow the optimum system transfer function in equation (8) to 
be written as a function of only Dj . Thus 




0.289d/'^s 


2.6d/‘^^s3 + 


+ 0.75Dj°‘^®s + 

3.4Dj°’^s2 + 2.6d/‘^^s + D 


( 11 ) 


where it will be recalled that 


Di = 




It is sometimes convenient to express equation (1 1 ) as 


V(s) 


0.289Di°-^ s2 + 2.6Di'^*OSs + 3.46Di°-l 


s'* + 2.6 Di°’^^s3 + 3.4d/’^s2 + 2.6d/*'^^s + 


( 12 ) 


Four observations can be made concerning this example. 

1 . The denominator is exactly of Butterworth form; that is, 

s'* + 2.6cos^ + 3.4o)^s^ + 2.6co^s + co'* (12a) 

This, in general, was anticipated because it has been shown that for systems with one input and a 
scalar Q, the QPI produces an optimum system having a Butterworth response (refs. 5 and 1 7). 
Equation (11) indicates that optimization of this particular system produces a Butterworth form 
for all values of q,, and r. It will be shown later that open-loop system dynamics greatly influence 
the magnitude of qj i/r required to produce a Butterworth response from the closed-loop system 
(see fig. 1 3). 



regulator cases for q/r = 10 and q/r - 1000. 


2. The transfer function resulting 
from the state tracking optimum systems 
always indicated unity dc gains for 
0A/0m see eq. (11)); however, 

transfer functions resulting from 
regulatory optimums did not produce 
unity dc gains (see ref. 5 and fig. 8). The 
feedback gain terms produced by the QPI 
for the regulatory optimums were the 
same as the corresponding terms for the 
state tracking, but the state tracking cases 
had one more feedforward gain term than 
the regulatory cases. This extra gain term 
is on the input and the input can be 
reformulated into other feedforward 
states plus acceleration. It was the 
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addition of the extra term on the attitude feedforward from the input that resulted in optimums 
with unity dc gains. Applications of the QPI requiring unity dc gain should not be optimized in the 
regulator form; however, the state tracking form can be used to obtain unity dc gain. 

3 . The numerator terms as expressed in the form given by equation ( 8 ) also yielded 
straight-line log-log plots. 

4 . The Butterworth numerator coefficient for N2 plotted on figure 7 is observed to be a 
good approximation for the computer derived values for N2 over this range of Dj . 

The calculations of this example will now be repeated but with a weighting factor on the rate 

error. 


Attitude and Rate Weighting Only 

When the system of figure 3 with f4i = f4 2 = f4 3 = f4 4 =0 was optimized for both attitude 
and rate weighting for several combinations of r, g4 1 , and with qi 1 = q2 2 0 -e., equal weights on 
both attitudes and rates), straight-line log-log plots for D4, D3, D2, versus Dj again resulted. 
Because of the effect of q2 2 on the system, the optimum coefficients for this case were no longer 
Butterworth. Nonetheless, as shown in figures 9 and 1 0 , the coefficients may be expressed by the 
same form, that is. 


D4 = 2 . 68 D ,°-2 8 

D3 = 3 . 55 Di«-”® 
D 2 =2.74Di°-8^ 



Figure 9.— Optimum numerator terms for 1/s^ plant 
(attitude and rate weighting). 


N, =D, 

N2=0.75Di‘-* 1 (13) 

N 3 =0.31Di‘-° ) 



plant (attitude and rate weighting). 
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Equations (13) provide the optimum follower transfer function for the attitude and rate 
weighting as follows: 


<J>m(s) 


O.SlDjS^ + 0.75Dj^‘^s + Dj 


+ 2.68d/’^®s3 + 3.55Dj°‘^^^s2 + 2.74Di°'®‘*s + 


The important feature here is that aU the terms of both preceding examples are characterized by a 
similar form of parametric relationship. All feedback and feedforward gains have a functional 
relationship to D, for attitude and rate weighting as well as attitude weighting only. 



Figure 11.— Attitude time response for I/s'* plant 
comparing attitude weighting to attitude and 
rate weighting. 


It should be noticed that the total 
outer-loop gain term, D, , is not influenced by 
the addition of Q 2 2 • However, as shown by the 
equations for Dj, D3, and D4 (eqs. (13)), 
q 2 2 does influence the inner-loop optimum 
coefficients and thus the dynamics of the 
follower response. This characteristic can 
perhaps be best shown by figure 1 1 , where the 
optimum system attitude time responses are 
shown for the follower transfer function of 
I/s'* that has been optimized by the QPI. In 
figure 1 1 , time histories of attitude and rate 
weighted cases (i.e., qj , ■= q 2 2 ) are visibly more 
damped than are the attitude weighted cases 
(i.e., q 2 2=0)- However, if qi 1 and q 2 2 are 
altered a faster responding system can be 
produced. Figure 1 1 indicates some of the 
flexibility within the QPI for producing various 
system responses. 


The calculations for the single-axis attitude weighting only example will now be repeated for 
the X— 14B aircraft system. 


X-14B System 

To extend the dependency of the feedforward and inner-loop feedback terms on Dj , consider 
a more physically realistic system in which not all the follower feedback terms f 4 i , f 4 2 > f 4 3 > and 
f 4 4 are zero. Such a system is shown in figure 3. The X— 14B example will be of the form of 
figure 3, but with f 4 , =0. 

The plant representation of the X— 14B, as described in Guidelines for Investigation and 
appendix B, consisted of an actuator and single-axis aircraft equation: 
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X 


<)>aCs) hjl 

^ _ 6 X 

u(s) ” s(s + Lp/Ix) 




4.585 

“ s(s + 0.45) "" 


62.8^ 

s2 + 2(0.5) (62. 8)s + 62. 8^ 


18,100 

s(s3 + 63.3 s 2 + 3978s + 1775) 


(14) 


where (L§/Ix)coa^ = g 4 1 is defined in appendix B. With the addition of system dynamics, the 
direct solution shown in the first example becomes considerably more laborious; thus only the 
graphical solution is presented here. The system transfer function will now be optimized for various 
attitude weightings. Since the system transfer function is fourth order, the optimum transfer 
function can be expressed in the form given by equation ( 8 ): 


^^(s) N3S2 + N2S + 

+ D 2 S + 

where, from table 1 , the total outer-loop gain is equivalent to Dj and 


Di =g4i(qi i/r)*/^ 

Di = 4.585(62.8)2 (qii/r)'/2 

Figure 12 indicates the manner in which the asymptotes for the D 2 , D 3 , and D 4 
denominator coefficients can be drawn. Equations (7e), (7f), and (7g) indicate that the exact values 
for D 2 , D 3 , and D 4 are; 


D2 “ §4 1 1^1 “1^4 2 
D3 “ 84 1 1^2 “ 1^4 3 

1^4 “ 84 1 1^3 “ 1^44 
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10 10 ^ 10 ^ 10 ^ 10 ^ 10 ^ 10 ^ 10 ® 10 ® 

D, — 


Figure 12.— Denominator asymptotes for X— 14B roll axis (attitude weighting). 

The asymptotes for the denominator coefficients D 4 , D 3 , and D 2 will for every small values of 
augmenting gains, be the open-loop system values of f 4 4 , f 4 3 , and Hence D 4 = 63,3, 
D 3 = 3978, and D 2 - 1775 are the asymptotes for this example. Next the asymptotes for 04 , 03 , 
and O 2 will be determined for the region where the augmenting gains are very large relative to the 
open-loop system values. Because the example uses attitude weighting only, as augmenting gains 
increase without limit, the optimal transfer function coefficients, O 2 , O 3 , and O 4 , approach the 
Butterworth coefficient asymptotes of 3.40i®'^, and 2.60i®'^^, respectively (see 

ref. 17). Between the small and the large augmenting gain regions, additional midrange asymptotes 
can sometimes be drawn. Midrange asymptotes can be drawn in our example because the 10— cps 
actuator dynamics are much greater than the other system dynamics and could be eliminated in the 
system representation for small augmenting gain values. Hence, cl>/^(s)/u(s) = + (Lp/Ix)l 

would be the open-loop transfer function being optimized if the 10 — cps actuator transfer function 
were considered to be unity. The asymptote for this optimized second-order open-loop transfer 
function for attitude weighting is a second-order Butterworth transfer function; that is, 




( 15 ) 


where, from table 1 , 


D/ = g4i(qi i/r)’/^ 
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To plot these second-ordei asymptotes for Dj on the fourth-order graph of figure 12 , 
equation (15) must be rewritten in terms of Di . Replacing the actuator of equation (14) by unity 
results in a numerator g 4 / of 

§4 1 “ 

hence 

g4 / = g4 1 /oJa^ 

and therefore the second-order D/ equals the fourth-order D, divided by , that is, 

D,' = Di/tUa^ 

Substituting this expression into equation (15) yields a transfer function in terms of Dj , 

4-^(5) Dj/W^2 

u(s) 2 , 2 -> 0-5 . 2 

s + 1.4(02/0)3'^) S + Dj/o)^ 


U)^s2 + 1.4 u)„d/*^s + D, 
a J- ^ 


3945s2 + + Dj 


Thus, the 88 D 1 term becomes the midrange asymptote for Dj. 

The preceding discussion indicates how asymptotes for denominator coefficients in the 
optimum transfer function can be drawn for aU values of Dj . The asymptotes in figure 12 should 
be a good approximation of the optimum coefficient values for attitude weighting only. Figure 13 
is a plot of the optimum solutions as well as the asymptotes previously shown in figure 1 2. It can be 
seen that over a considerable range of the variable Dj , the asymptotes provide a good 
approximation of the true optimums. 
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Lg/I, = 4.585 
Lp/ 


Figure 13.- Denominator terms for X-14B roll axis (attitude weighting). 


Bulier worth 


2nd order 




Figure 14.— Numerator terms for X— 14B roll axis (attitude 
weighting). 

An asymptote for the optimum numerator coefficient N 2 can be found in much the same 
way as the denominator coefficients were determined. Figure 14 is a plot of the numerator 
coefficient asymptotes and the actual optimum values. 









CONCLUSIONS 


The quadratic performance index was found to be the most versatile criterion for the design 
of a prefilter model reference system for a VTOL aircraft. Calculation of the optimum system gains 
for this criterion can be determined through readily available digital computer programs. The 
system, which was optimized in this study, consists of a second-order prefilter model and a 
simplified fourth order (single input, single control, constant numerator, and non- type- zero) roll 
axis representation of the X— 14B aircraft. 

Expressions for the optimum total outer-loop gain were analytically determined for both type 
zero and non -type -zero systems. Analysis of the canonical form of the state equation and the 
Riccati equation for single-input, single-control systems indicates that the total-outer-loop gain can 
always be analytically determined. For type-zero systems the total- outer-loop gain includes the 
outer -loop weighting element of the Q matrix (the 1 — 1 element in the report examples), the 
control weighting element of the R matrix, the total system gain element, and the outer-loop gain 
term of the unaugmented follower. As the outer- loop gain term of the unaugmented follower 
approaches zero the optimum total-outer 4oop gain approaches that of the non-type-zero system, as 
was expected. For the non-type-zero system the total -outer -loop gain is the product of the total 
system gain and the square root of the quotient of the 1 — 1 element of the Q matrix divided by the 
scalar element of R (only single-input, single-control systems were studied). 

The only terms that affect the optimum outer-loop feedback gains are the attitude weighting 
term, the control weighting term, the total system gain term, and, in the case of the type-zero 
system, the outer- loop gain term in the unaugmented follower. Hence the addition of other 
weighting terms (such as rate weighting) does not affect the outer-loop feedback gain. Weighting 
terms other than on attitude will, however, influence the inner-loop dynamics of the follower. 

A graphical method for approximating the optimal gains was developed. By this method 
asymptotes were generated for the inner -loop feedback gains in terms of the total outer-loop gain 
over a wide range of weighting. The asymptotes approximate the loop gain values for the inner 
feedback and the feedforward terms. For the X— 14B example, with attitude weighting only, the 
asymptotes plotted on log-log coordinates are piecewise linear. For small weighting of attitude the 
asymptotes represent the characteristics of the unaugmented X— 14B, and for large weighting of 
attitude the asymptotes are of the fourth-order Butterworth form. For this particular example, an 
additional set of mid-range asymptotes result which are of the second-order Butterworth form. With 
the addition of equal rate weighting in the performance index, straight-line asymptotes can again be 
constructed - however, not of Butterworth form. The asymptotes for this and other combinations 
of weighting are constructed with the aid of a digital computer. 
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Applications of the quadratic performance index requiring unity dc gain should not be 
optimized in the regulator form; however, the state tracking form can be used to obtain unity dc 
gain. 

Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif. 94035, Aug. 18, 1970 
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appendix a 


HISTORICAL REVIEW 


This appendix briefly reviews the background of three areas of research investigated in this 
study: (1) characteristics of different performance criteria, (2) desirable handling qualities models, 
and ( 3 ) model reference systems. 

Various performance criteria were studied in references 8 to 1 3 . In principle these criteria can 
be used to optimize any order or type system; however, their effectiveness depends on whether or 
not the chosen error criterion is the one to be minimized. What is considered optimum in one 
application mi^t not be considered optimum in another. This section of the appendix will discuss 
some of the considerations for using performance criteria for model reference systems rather than 
for the conventional response feedback system. 

A prefilter model follower system may, in many ways, be considered an extension of a 
conventional response feedback system. To illustrate this extension, consider a single-input, 
single-output response feedback system. Such a configuration can be considered a model follower 
with a unity model, in which the system (or follower) is attempting to follow the model (the 
input). In other words, the response feedback configuration contains a model with infinite 
bandwidth and thus a follower to model bandwidth ratio of zero. A model with a less than infinite 
bandwidth would then introduce a nonunity transfer function as a model; hence the scheme would 
take on a more conventional model following form. A decrease in the model bandwidth would 
increase the follower to model bandwidth ratio and hence improve the ability of the follower to 
match the model’s response. 

Because of this rough similarity, performance criteria often cited in response feedback studies 
were the first to be considered for application to model reference systems. The major performance 
criteria discussed in these studies are: 


Ii = J edt 

Q 

Integral of error 

I2=J ie|dt 
0 

Integral absolute error 

rT 

I3 = J dt 

n 

Integral error squared 

II 

?fr 

0. 

r-b 

Integral time absolute error 

0 

Is = J tle^ dt 

Integral time squared error squared 

0 

Is=J F[e,t]dt 
•^0 

Integral weighted error 
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Generalized performance index 


h 



+ tj + . . . + t^”(d%/dt*^)^Jdt 


Ig = J^le^Qe + u^Ruldt 
0 


Quadratic performance index 


Each criterion had some advantages and some deficiencies, and the proper criterion for a given 
system optimization should be chosen in light of its characteristics. Criteria I, through Ig are 
usually applied to single error terms. Criteria I2, I4, and some forms of Ig and I7 are difficult to 
evaluate analytically. Criterion Ig was chosen because it includes the cost of control, produces a 
closed form solution for a multiparameter optimization integral, and can be expeditiously 
optimized through a digital computer program. The indicial response characteristics for response 
feedback systems for each of these criteria is discussed in the following paragraphs. For criteria I2 , 
I3, I4, and Ig comments regarding optimums produced for model reference systems are included. 
The prefilter model is that described in the section on desirable handling qualities. 

Criterion is quite adversely affected when used to optimize systems that overshoot 
because the overshoots decrease rather than increase the integral. This criterion cannot be used 
satisfactorily in producing an optimum second-order system. (The minimum Ij is for zeta = 0 . 0 .) 

Criterion I2 has moderate selectivity in optimizing systems up to sixth order, and is intuitively 
appealing since it represents the area of error between the model and the follower time histories. 
Optimization through the use of this criterion leads to a minimization of this area. While absolute 
value functions are very difficult to use analytically, they are easily mechanized on analog or digital 
computers. For the second-order system, the optimum f is 0 . 7 . When used with model reference 
systems, I2 produced very satisfactory, nonoscillatory follower responses. Followers of second, 
fourth, and sixth order were studied with repetitive operation analog computers. 

Criterion I3 is the most widely used of aU the criteria primarily because of its mathematical 
convenience; however, systems thus optimized tend to be more underdamped than do systems 
optimized by other criteria. I3 can be used with statistical inputs and Parseval’s integral technique 
has been used in evaluating this integral exactly. Variations of this integral, such as Is, can be 
evaluated by Parseval’s technique; however, system orders greater than five are very cumbersome to 
evaluate. When used with model reference systems I3 produced unsatisfactory follower responses. 
The optimum I3 gains produced follower responses that exhibited the same oscillatory nature as 
were exhibited for the response feedback systems. Hence this criterion was rejected for application 
on the X— 14 B aircraft. 

Criterion I4 has been lauded by Graham and Lathrop (ref. 8) as a criterion with better 
selectivity than I2 or I3 , and one that picks “good” systems; however, it is difficult to evaluate 
analytically. For a response feedback system I4 has been found most suitable when optimizing to 
a step input (see ref. 8). The time weighting in the integrand tends to reduce the cost of the large 
initial error between the step input and the follower output, which is inherent in the response 
feedback system. However, for a model following system no such large initial error should exist 
between the model and the follower. In fact, for a step input into the model, the follower should 
follow the model “closely” at all times. Thus, for such a system, the time weight feature 
of I4 serves no useful function and hence its use was not considered further. 
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Criterion 1 5 yields step responses that are less oscillatory than I3, and compares favorably 
with the more desirable system optimums provided by I4. While I5 is more cumbersome than I3, 
it can still be evaluated analytically. 

Criterion Ig is a comprehensive representation for any of the above criteria and has extra 
provisions for special error weightings to be chosen according to the particular design objectives. 
For example, special functions might penalize overshoots more than undershoots. 

Criterion I7, a generalized performance index proposed by Aizerman (ref. 18), contains the 
system error, its first n time derivatives, and the constants tj, where i = 1 , 2, . . . n, determined 
from the differential equation of the desired system response. This criterion is included because it 
leads into Ig and gives some insight into the choice of the matrix elements of Ig. 

Criterion Ig is the quadratic performance index explained in appendix B. It has several 
significant advantages, one being that it includes the cost of control in the integral. The ability to 
weight the various states of the solution enables the designer to choose the proper weighting 
matrices to provide the particular desired responses. That is, if the obtained response is too 
oscillatory, then through the weighting matrix element changes the desired, less oscillatory, 
response can be obtained. The quantity Q/R should be chosen according to the desired bandwidth 
ratio of the follower to the model. The disadvantages of the QPI include the following; the index 
must be expressed only in terms of states and the control vector, all states must be observable, only 
linear system formulations arc acceptable. Even with these disadvantages, it is felt that Ig is the 
best performance criterion. 

Various desirable models were studied in references 1 , 3, 6, 7, 1 9, and 20. The model chosen 
for this study was the roll axis for a VTOL aircraft which had been optimized for handling qualities 
in hover. The model, expressed in terms of the desired aircraft roll attitude, represents the optimum 
attitude control system determined on a large moving base simulator for a variety of tasks (see 
ref. 7). The criterion optimized in reference 7 to obtain the desired model was pilot opinion; hence 
an implicit and unknown weighting was assigned by the pilot to roll attitude, rate, acceleration, 
lateral translation, stick forcing functions required, controlled element dynamics, and the 
manipulator. Other variables that affected the rating were the visual scene, g-levels, training, fatigue, 
motivation, illumination, vibration and temperatures. In reference 7, the optimum attitude system 
in roll was found to require a damping of —2.8/sec and a natural frequency of 2 rad/sec. As a 
matter of interest, this system also required a control power in excess of 1.2 rad/sec^ and a stick 
sensitivity of 0.5 rad/sec^ /in. This system yielded an average pilot rating of 2 on the Cooper Pilot 
Rating Scale (see ref. 21 ). Simulated conditions were for calm air; hence to extend results to flight, 
additional control power must be included for trim and upset disturbances. The roll manipulator 
had a maximum control deflection of ±5 inches with a force gradient of 1.8 Ib/in. and a breakout 
force of 1 lb. 

For comparison, the hover model used for the Bell X— 22 aircraft was co = 3 rad/sec, 
and f = 1.0 (ref. 6). The pitch axis in the X— 15 uses a model of co = 2 rad/sec and f = 0.7 for the 
rate command (ref. 1 9). 

Model reference systems are discussed in references 1, 3, 4, 5, 6, 22, 23, and 24. Several types 
of model reference systems have been discussed in the literature, primarily the 
following: (1) prefilter model following systems (sometimes called the model in the system or the 
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high gain system), (2) implicit model or model in the performance index systems, (3) inverse model 
systems, and (4) parameter adjustment systems. All of these approaches with the exception of (2) 
are considered to be adaptive, with the gains either slewed or fixed. The intended application of a 
model reference system is the chief determinant of the best system configuration. The following 
paragraphs review some characteristics of four types of model reference systems. 

The configuration under consideration in this study is the prefilter model following system 
with fixed gains. The follower used in these examples was a VTOL hovering aircraft. It was also 
assumed that an onboard airborne computer on which to program the model would be available. 

The prefilter model following system requires a physical or programmed model to which the 
controlled system or follower alines itself. If the system to be optimized is to operate over a wide 
range of conditions, then it may be desirable for the follower to be altered. An example of a 
prefilter model following system using this approach is the NASA-FRC Lockheed Jetstar (refs. 1, 
3, 4, 22) in which dynamic pressure can be used to vary the appropriate gain terms on the errors 
between the model and corresponding aircraft parameters. Varying or slewing of gain need not be 
restricted to the follower but may include a variable model as well. A possible application of this 
most general case is a VTOL aircraft going from hover to cruise, in which case the desirable model 
to be simulated is continuously changing. 

One primary limitation of the prefilter model following system is the operational problem 
associated with obtaining sufficient closed loop follower bandwidth. Factors affecting bandwidth 
include sensor and system noise, structural vibration feedback, nonlinearities, and sensor and 
actuator dynamics. These effects, however, will also degrade the performance of the other 
approaches. 

For “good” following, the follower bandwidth should be roughly twice the model bandwidth 
(refs. 1 , 22). However, the high-frequency gust response effects on the pilot for this type of system 
are not clear. Although these effects have not been a problem with this type system in the X— 15 
(ref. 19), it should be noted that the X— 15 disturbance content is not as high as that for a VTOL. 
Another limitation of the model following system is that the complexity of the system to be 
modeled generally requires the use of an airborne computer. 

An advantage of the prefilter model following system is its ability to achieve both static and 
dynamic matching of the follower to the model. Other desirable features of this system are its 
ability to null the effects of a follower with either poorly defined or slightly variable parameters 
and its inherent ability to isolate or decouple the axis of control. 

The model in the performance index is a nonadaptive technique that provides a way of 
matching desired dynamics. In implementing this technique, the model is not available for reference 
when the system is operating, because the system has already been optimized and the model 
removed. For example, if changes in the plant occur after the system has been optimized, then the 
uncorrected effect is felt on the plant outputs. The advantages of this technique are simplicity and 
low cost of implementation. For instance, a computer on which to program the model is not 
required. 

The inverse model in the feedback approach uses a high gain and some compensation in the 
forward loop so that the closed-loop transfer function approaches the ratio of the forward transfer 
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function over the open-loop transfer function. In the limit, the closed-loop transfer function 
becomes the reciprocal of the feedback transfer function (which is the inverse of the model). Thus 
the closed-loop system becomes the desired model. The model in this approach usually has more 
poles than zeros; hence the inverse requires at least one differentiation that introduces noise into 
the system. An operational flight test of a high gain saturating version of this approach at Ames on 
an F— 102 was successful (ref. 25). The test was conducted for Mach numbers from 0.36 to 1.15 
and altitudes from 1 0,000 to 40,000 feet. The system limited longitudinal natural frequencies to 
one-half the basic aircraft frequencies, and the damping ratios were improved to a range from 0.56 
to 0.69. However, objectionable features were sensitivity to gust disturbances and high 6 to 
5 elevator gains at low speeds. 

In the model reference adaptive parameter adjustment technique, an adjustment mechanism is 
used to alter compensation elements so that the form of the overall system is the same as that of 
the model. The parameter adjustment techniques can be grouped into three general categories in 
which adjustment of the compensation systems is based on the gradient of the error functions, 
Lyapunov’s direct method, and partial derivatives of functions. In all these approaches, a nonzero 
input is necessary for adjustment. With no inputs, errors between the actual and desired values can 
result. In both the first and third categories, the error function has a large effect on the proper 
operation of the adaptive portion of the system. 
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APPENDIX B 


MATHEMATICAL FORMULATION 


The techniques of state space modeling and of minimizing the quadratic performance index 
are detailed within this section. 


State Space Formulation 

Figure 1 5 illustrates a composite block diagram of the type system considered in this report. 
The follower is expanded into three cascaded second-order transfer functions (i.e., actuator, 
bending mode, and linearized roll hover equation). In this diagram, each gain block (e.g., Kjjj, RC’ 
Kr, represents an array of scalar gains that operates on aU of the states of the model and 
follower. When devoid of the bending mode, this system reduces to the fourth-order configuration 
considered in the test. A further simphfication such as approximating the actuator by a unity gain 
would reduce the follower to second order and thus allow for the application of perfect model 
following techniques (refs. 1 5 and 24). 



Figure 15.- Sixth-order follower block diagram. 


There are several possible mathematical formulations of the system shown in figure 15. Hence 
the development of the particular representation utilized in the text for the fourth-order example 
cited above is given below. Consider this fourth-order follower system as configured in figure 1 6. 
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Follower 



Figure 16 — Detailed system diagram. 


The corresponding system equations are; 


Model Equation 


(Bl) 

Follower Equations 

coa^u = (s^ + 2faWaS + <^a^)^x 

(B2) 


k(L5/Ix)5x = s[s + (Lp/Ix)]^A 

(B3) 

Control Equation 

u = -(kgj^s + kg^)6x - + k^^)0A + <t>c 

(B4) 

where 

^c = (k0mS + k0^)0ni + k55 


An equivalent configuration representing the follower in its canonical form 
first combining equations (B3) and (B4), that is, 

can be generated by 





Then combining equations (B2) and (B3) yields a composite follower equation from which the 
equivalent or representative system can be sketched: 



= - ^4453 - f^3s2 - f^2^)4>A (B 6 ) 

from which evolves figure 1 7. 



Figure 17.- Representative system diagram. 


This figure represents a special case (i.e., f 4 i =0) of the general representation shown on figure 3 
in the text. The corresponding equations for figure 1 7 are 
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Model Equation 




Follower Equation 

g41 U = [S'^ - f44S^ - f43S^ ~ f42S]0A 


where 


§4 1 ~ k( Lg/Ix)i^a^ (B7) 

Control Equation 

u = [-kjS^ - kjS^ — k,s - ko]0A + <^c 

From these equations and again for 5 equal to a step input, the state equation is 
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where 

§4. = k(L5/Ix)coa^ 

It is worth noting that equation (B9) is a special example of equation (5) in the text. For this 
representation, the total outer-loop gain Dj referred to in the text is 

D. =g4i(q.i/r)‘/" =k(L5/Ix)oja^(q,,/r)‘/2 

The nonhomogeneous state equation (eq. (B9)) can be programmed directly onto the analog 
computer. However, for digital computation, it is highly desirable to operate with a set of 
homogeneous equations. Under certain conditions (when the input can be represented by a state 
equation), such a set of equations can be generated from the nonhomogeneous form; a discussion of 
this procedure foUows. 

Homogeneous State Equation 


X = FX Xo = X(0) 


(BIO) 
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For systems in which the control equation can be expressed entirely in terms of the system states, 
that is, 

u=KX 


the homogeneous state equation can be written as 


d 

dt 


X' 

6 


X 



+ [GK] 


(BID 


The model following systems considered in this study comply with this condition. With the 
fourth-order follower example shown in figure 18 for illustration, the homogeneous state equation 
of the form given by equation (BIO) becomes 
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0 

0 _ 

_ 6_ 


The homogeneous equation for the second- and sixth-order examples may be generated in a similar 
manner and will not be shown here. 

There are two factors worth noting in the preceding development. First, the natural or basic 
outer-loop follower feedback term is equal to zero (i.e., f 4 1 = 0). Thus the transfer function for this 
follower is categorized as a non-type zero. The second factor is the distinction between the system 
matrix F of equation (B9) and F of equation (B1 !)• F consists of the basic unaugmented 
follower and model matrices Ff and M, divorced from the augmenting feedback GK; that is. 


Ff 

0 1 

0 

0 

M 1 

B 

0 

0 1 

E 


(B13) 
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On the other hand, F represents a composite system matrix incorporating the augmenting feedback 
gains GK; that is, _ 

F = F + GK (B14) 

Hence, in the text, F and Ff are referred to as the unaugmented system matrix and the 
unaugmented follower matrix, respectively, and F is referred to as the augmented system matrix. 

The goal achieved by the preceding development was the illustration of the modeling 
techniques used in this study and further, for ease of solution, the formation of a homogeneous 
state equation from a nonhomogeneous form. 


Quadratic Criterion 

The quadratic performance criterion, as given below, forms the basis of much modern and 
optimal control theory. Its minimization is achieved through a variety of techniques, including 
Hamilton-Jacobi, calculus of variations, or Pontryagin’s maximum principle, most of which evolved 
into the solution of a first-order matrix differential equation referred to as the matrix Riccati 
equation. For a concise development of the Riccati equation, the reader is referred to references 5, 
6, and 9. The Riccati equation, whose solution provides the designer with a set of system gains 
that will minimize the quadratic criterion, is given as 

P = PGR-‘ G^P - PF - F^P - hTqh (B1 5) 

A formulation and usage of this optimal technique will now be developed. 

The model reference system under investigation is the prefilter type. To delineate the notation 
for the optimization of the regulator problem (i.e., no driving functions, 6 = 0), consider again the 
second-order model and linearized single-axis second-order aircraft follower configuration: 



Control Follower (Aircroft) 


Figure 1 8 .— Second-order follower system. 
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Statement of problem— Find a control u through a selection of gains K (i.e., ^0m’ 

k^a> ^ 0 ^) such that the quadratic performance index is minimi 2 ed. 


V = / xTrTqHX + uTru dt QPI 
0 


(B16) 


Solution to the problem— Given the matrices F, G, H, Q, and R, which will be defined 
subsequently, the specific gains, K*, which minimize the QPI, can be determined through the 
solution of the matrix Riccati equation by the expression 

K* = -R-igTp (B17) 

where P is the solution to the Riccati equation (eq. (B15)). The optimum control becomes 

U* = — R‘‘ g"^PX Optimum Control 

To define the notation used in this optimization procedure, consider the modeling of the 
second-order follower system shown above. 
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— 1 
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dt 
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Model Equation 
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0A 




Control Equation 


u = KX 
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Follower Equation 
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Xf = Ff X£ 
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The key step in the optimization of this system is the determination of a matrix of gains K, 
which minimize the QPI. These gains can be determined through equation (B17), which requires 
first the solution of the Riccati equation in order to establish P. The Riccati equation is rather 
unruly; however, it can be handled nicely on a digital computer utilizing existing subroutines of the 
Automatic Synthesis Program (ASP) or by a more recently developed Fortran IV version of ASP, 
referred to as FASP (ref. 16). Required as inputs to these subroutines are the matrices F, G, H, Q, 
and R, which for the examples of figure 1 8 are shown below. The system matrix F is formulated 
as 



System Equation 


where the gains K in our case were initialized at zero. The other matrices are 


X = 




'm 


G = 


Gf 

0 


<I>A 
<^A 
‘I’ln 

V 

0 

Mix 

0 

0 


H = 


0 10 


0-101 


Q = 


^11 

0 


0 ■ 
^22 


R = r 


11 


Matrices H, Q, and R, all of which appear in the QPI integrand, are arbitrary matrices selected 
by the designer; H is constructed so as to include in the integrand the errors which are to be 
minimized. In the case above, the differences between the model and aircraft attitudes and rates are 
to be minimized; that is. 


HX 





‘t’A 

^A 


Matrices Q and R are the selected weighting matrices on the errors and on the cost of control. For 
this example, the scalar qj j is the weighting on the error of the attitude match, qjj on the rate 
match, and ri j on the control u. 
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This optimization technique, by its nature, requires that all states be available for feedforward 
or feedback. Such is, of course, not always the case; however, in many instances a near optimum 
can be achieved when a partial set of states is available. An additional feature of this procedure is 
that a different set of gains result from each combination of H, Q, and R. Hence there is a need to 
better understand the effects of changing these matrices on the optimums they produce. 

Once the above system matrices have been generated as inputs for the programs ASP or FASP, 
the steady-state Riccati matrix P and the optimal feedback gain matrix K* can be determined by 
equation (B17). This optimal gain matrix can then be used to construct a model reference system 
that is “optimized” for the imposed set of constraints. The P matrix from which the K matrix is 
derived also contains much information regarding the ability of the optimized system to regulate 
about a variety of initial conditions or upsets. This information evolves from the expression 


where 


^min ~ ^o^^o 


Xo=X(0) 


(B 18 ) 
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